Casimir densities for two concentric spherical shells 
in the global monopole spacetime 



A. A. Saharian ^*and M. R. Setare ^ ^ 
^ Department of Physics, Yerevan State University, Yerevan, Armenia 

and 

^ Institute for Theoretical Physics and Mathematics, Tehran, Iran 
Department of Science, Physics Group, Kurdistan University, Sanandeg, Iran 
Department of Physics, Sharif University of Technology, Tehran, Iran 

February 1, 2008 



Abstract 

The quantum vacuum effects are investigated for a massive scalar field with general 
curvature coupling and obeying the Robin boundary conditions given on two concentric 
spherical shells with radii a and b in the D + 1-dimensional global monopole background. 
The expressions are derived for the Wightman function, the vacuum expectation values of 
the field square, the vacuum energy density, radial and azimuthal stress components in the 
region between the shells. A regularization procedure is carried out by making use of the 
generalized Abel-Plana formula for the series over zeros of combinations of the cylinder 
functions. This formula allows us to extract from the vacuum expectation values the parts 
due to a single sphere on background of the global monopole gravitational field, and to present 
the "interference" parts in terms of exponentially convergent integrals, useful, in particular, 
for numerical evaluations. The vacuum forces acting on the boundaries are presented as 
a sum of the self-action and interaction terms. The first one contains well known surface 
divergences and needs a further regularization. The interaction forces between the spheres 
are finite for all values a < b and are attractive for a Dirichlet scalar. The asymptotic 
behavior of the vacuum densities is investigated (i) in the limits a — > and & — > cxd, (ii) in 
the limit a, 6 — > oo for fixed value b — a, and (iii) for small values of the parameter associated 
with the solid angle deficit in global monopole geometry. We show that in the case (ii) the 
results for two parallel Robin plates on the Minkowski bulk are rederived to the leading 
order. 
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1 Introduction 



The Casimir effect is regarded as one of the most striking manifestation of vacuum fluctuations 
in quantum field theory. The presence of reflecting boundaries alters the zero-point modes of a 
quantized field, and results in the shifts in the vacuum expectation values of quantities quadratic 
in the field, such as the energy density and stresses. In particular, vacuum forces arise acting 
on constraining boundaries. The particular features of these forces depend on the nature of the 
quantum field, the type of spacctimc manifold and its dimensionality, the boundary geometries 
and the specific boundary conditions imposed on the field. Since the original work by Casimir 
in 1948 [1] many theoretical and experimental works have been done on this problem (see, 
e.g., [2, 3, 4, 5, 6, 7, 8, 9] and references therein). Many different approaches have been used: 
mode summation method with combination of the zeta function rcgularization technique, Green 
function formalism, multiple scattering expansions, heat-kernel scries, etc. The Casimir effect 
can be viewed as a polarization of vacuum by boundary conditions. Another type of vacuum 
polarization arises in the case of external gravitational fields [10, 11]. In a previous paper [12] 
we have studied an example of a situation when both types of sources for the polarization are 
present. Namely, we have investigated the vacuum expectation values of the square of a scalar 
field and energy-momentum tensor induced by a spherical shell in the spacetime of a point-like 
global monopole. 

Topological defects have attracted a great deal of attention because of their relevance to a 
number of different areas ranging from condensed matter to structure formation (for a review 
see [13]). In the context of hot big bang cosmology, the unified theories of the fundamental 
interactions predict that the universe passes through a sequence of phase transitions. These 
phase transitions might have given rise to several kinds of topological defects depending on the 
nature of the symmetry that is broken [14]. If a global SO (3) symmetry of a triplet scalar field is 
broken, the point like defects called global monopoles are believed to be formed. The simplified 
global monopole was introduced by Sokolov and Starobinsky [15]. The gravitational effects of 
the global monopole were studied in Ref. [16], where a solution is presented which describes a 
global monopole at large radial distances. The quantum vacuum effects of the matter fields on 
the global monopole background have been considered in [17, 18, 19, 20]. The effects produced 
by the non-zero temperature are investigated as well [21]. The zeta function and the Casimir 
energy for spherical boundaries in this background are considered in [22, 23]. 

In this paper, we continue the investigation started in [12] and consider the Casimir densities 
in the region between two concentric spherical shells on background of the D + 1-dimensional 
spacetime of a point-like global monopole. The positive frequency Wightman function, the 
vacuum expectation values of the field square and energy-momentum tensor are investigated for 
a massive scalar field with general curvature coupling parameter ^. In addition to describing the 
physical structure of the quantum field at a given point, the energy-momentum tensor acts as 
the source of gravity in the Einstein equations. It therefore plays an important role in modeling 
a self-consistent dynamics involving the gravitational field [10]. We study the general case of 
Robin boundary conditions with different coefficients for the inner and outer spheres. The Robin 
boundary conditions are an extension of the ones imposed on perfectly conducting boundaries 
and may, in some geometries, be useful for depicting the finite penetration of the field into 
the boundary with the "skin-depth" parameter related to the Robin coefficient [24, 25]. It is 
interesting to note that the quantum scalar field satisfying Robin condition on the boundary 
of a cavity violates the Bekenstein's entropy-to-encrgy bound near certain points in the space 
of the parameter defining the boundary condition [26]. This type of conditions also appear in 
considerations of the vacuum effects for a confined charged scalar field in external fields [27] 
and in quantum gravity [28, 29, 30]. Mixed boundary conditions naturally arise for scalar and 
fermion bulk fields in the Randall-Sundrum model [31]. In this model the bulk geometry is a 
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slice of anti-de Sitter space and the corresponding Robin coefficient is related to the curvature 
scale of this space. For scalars with general curvature coupling the essential point is the relation 
between the mode sum energy, evaluated as a renormalized sum of the zero-point energies for 
each normal mode of frequency, and the volume integral of the renormalized energy density. In 
[32] it has been shown that in the discussion of this question for the Robin parallel plates it is 
necessary to include in the energy a surface term concentrated on the boundary (see [33, 34, 35] 
for similar issues in more complicated geometries and the discussion of the paper [32] in [36]). 

We have organized the paper as follows. In the next section we derive a formula for the 
Wightman function in the region between two spheres. The reason for our choice of the Wight- 
man function is that this function also determines the response of the particle detectors in a 
given state of motion. Following Refs. [33, 34, 35, 37], to evaluate the bilinear field products 
we use the mode sum method in combination with the summation formulae from [38] (see also 
[39]). These formulae allow (i) to extract from vacuum expectation values the parts due to a 
single sphere on the global monopole background, and (ii) to present the "interference" parts 
in terms of exponentially convergent integrals involving the modified Bessel functions. The vac- 
uum expectation value of the field square is obtained computing the Wightman function in the 
coincidence limit, and is investigated in section 3. The various limiting cases arc considered in- 
cluding the limit of the strong gravitational field corresponding to small values of the parameter 
associated with the solid angle deficit in the global monopole spacetime. Section 4 is devoted 
to the vacuum expectation values of the energy density and stresses. The formulae are derived 
for the interaction forces between the spheres. The behavior of these quantities is investigated 
in various limiting cases. Section 5 concludes the main results of the paper. 



2 Wightman function 

In this paper we consider a real scalar field (p with curvature coupling parameter ^ in the D + 1- 
dimensional spacetime of a point-like global monopole. In the hypcrsphcrical polar coordinates 
(r, (p) = (r, 61,62, ■■ ■ 6n, (p), n = D — 2, the corresponding geometry is described by the line 
element 

ds'^ = df - dr^ - a^r'^dnl, (1) 

where dftj^ is the line element on the surface of the unit sphere in D-dimensional Euclidean 
space, the parameter a is smaller than unity and is related to the symmetry breaking energy 
scale in the theory. In the spacetime given by line element (1) the solid angle deficit is {l — a'^)So, 
with Sd = 27r^/Vr(I?/2) being the total area of the surface of the unit sphere in D-dimensional 
Euclidean space. The field equation has the form 

(ViV^ + m^ + ^R)<p = 0, (2) 

where m is the mass for the field quanta, Vj is the covariant derivative operator associated with 
the metric given by line element (1), and 

R = n{n + 1)^^^^ (3) 

is the corresponding Ricci scalar (we adopt the convention of Birrell and Davies [10] for the 
curvature tensor). Note that for a 1 the geometry is singular at the origin (point-like 
monopole), r = 0. In (2) the values of the curvature coupling parameter ^ = 0, and 4 = with 

= {D — 1)/4Z) correspond to the minimal and conformal couplings, respectively. 

In this paper we are interested in the vacuum expectation values (VEVs) of the field bilinear 
products on background of the geometry described by (1), assuming that the field satisfies the 
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Robin boundary conditions 



- d 

Ar + Br-^ ] ^{x) = 0, r = a, fe, (4) 



on two spheres with radii a and b, a < b, concentric with the monopole. Here the coefficients Aj. 
and Br are constants, in general, different for the inner and outer spheres. The imposition of 
this boundary condition on the quantum field ip{x) leads to the modification of the spectrum for 
the zero-point fluctuations and results in the shift in VEVs for physical quantities. In particular, 
vacuum forces arise acting on constraining boundary. This is the familiar Casimir effect. The 
VEVs for the physical quantities bilinear in the field can be evaluated if the corresponding 
Wightman function is known. For this reason we first concentrate on the positive frequency 
Wightman function G+(x, x') = {0\ip{x)ip{x')\0), where |0) is the amplitude for the corresponding 
vacuum state. The Wightman function also determines the response of the particle detectors in 
a given state of motion [10] . By expanding the field operator over eigenfunctions and using the 
commutation rules one can see that 

(0|<^(x)^(a;')|0) = M^)^^^^'), (5) 

a 

where {ipa{x), ^a{x')} is a complete orthonormal set of positive and negative frequency solutions 
to the field equation with quantum numbers a, satisfying boundary conditions (4). Note that 
for D = 1 we have the standard Casimir-like geometry on Minkowski spacetime. 

In the hypcrspherical coordinates, for the region between two spheres the complete set of 
solutions to (2) with scalar curvature from (3), has the form 

iPa{x) = Par~''/''gui{Xa,\r)Y{mk;^,(f>)e-'^\ Z = 0,1,2,..., (6) 
where = (mo = I, mi, . . . , m„), and mi, m2, . . . , m„ are integers such that 

< nin-i < m„_2 < • • • < mi < I, -m^-i < nin < nin-i, (7) 

Y(mi;;i},(p) is the surface harmonic of degree I (see [40]). In (6) 

g,, (Aa, Ar) = J,, i\r)Y(^^ (Aa) - J^f {Xa)Y,, (Ar), A = VuJ^-m\ (8) 

J,y{z) and Yi,{z) are the Bessel and Neumann functions, and the functions with overbars are 
defined in accordance with 

f(°)(z) = AaF{z) + BazF'{z), Aa = Ac,- B^n/2a, B^ = B^/a, a = a,b. (9) 

Here and below we use the following notation 

1 



1^1 = - 
a 



1/2 

(10) 



assuming that is non-negative. This corresponds to the restriction on the values of the 
curvature coupling parameter for n > 0, given by the condition 

2 

na 

^- ~4(n + l)(l-(T2)- ^ ' 

This condition is satisfied by the most important special cases of the minimal and conformal 
couplings. The coefficients Pa in (6) can be found from the normalization condition 



/ 



Mx)fV^dV=^, (12) 
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where the integration goes over the region between the spheres, a < r < b. Substituting 
eigenfunctions (6), and using the relation 



/ 



\Y{mk;^,(t>)rdn = N{mk) (13) 



(the exphcit form for N{mk) is given in [40] and will not be necessary for the following con- 
siderations in this paper) for the spherical harmonics and the value for the standard integral 
involving the square of a cylinder function [41], one finds 



^2 _ vr^AT^f (Va,Ao) 



Pa = .... ' 'n- > (14) 



where we use the notation 



rj = -. (15) 



The functions chosen in the form (8) satisfy the boundary condition on the sphere r = a. Prom 
the boundary condition on r = 6 one obtains that the corresponding eigenmodes are solutions 
to the equation 

C^f (6/a, Aa) = 4«) (Aa)yj^) (A6) - (A6)yJ") (Aa) = 0. (16) 

Below the roots to this equation will be denoted by ju^k = Aa , k = 1,2,.... The corresponding 
eigenfrequencies ui = ijJvi,k are related to these zeros as tjJvi^k = \J A^^ ^./a^ + rv?. 

Substituting the eigenfunctions into the mode sum (5) and using the addition formula 

E T77^^(mfe;^,^)y*(m,;^^0') = '^l±I^C'^l\cosd), (17) 
for the expectation value of the field product one finds 

2( l\-n/2 ~ 

{()\^{x)^{x'm = ^A_L_ ^(2/ + n)Cr/'(cose)^M7.,,fc)r;,^(Va,7..fc), (18) 

^ 1=0 k=l 

with the function 

h{z) = ^=^==5.,(^, W«)5.,(^,^r7a)e*^/^V^^(^'-^), (19) 



and Cf{x) is the Gegenbauer or ultraspherical polynomial of degree I and order p. To sum over 
k we will use the generalized Abel-Plana summation formula [38, 39, 33] 



2 00 



h{x)dx 



Ki''^ (tjx) [h{xe^'/^) + /i(xe-^^/2)j 



Ri''^ {x) Ri"^ ix)li'^ (rix) - k'^^ irjx)li''^ (x) ' 

where Iu{x) and Ki,{x) arc the Bessel modified functions. The corresponding conditions for the 
function (19) are satisfied if r + r' + — < 26. Note that this is the case in the coincidence limit 
for the region under consideration. Note that we have assumed values Aa and Ba for which 
all zeros for (16) are real. In case of existence of purely imaginary zeros we have to include 
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additional residue terms on the left of this formula (see [38, 39, 33]). Applying to the sum over 
k in (18) formula (20), for the corresponding Wightman function one obtains 



- r K^^hn^) cy(.,../a)d°)(.,.//a) ^^^^ \^^rj^.^,> _ ^ 

where we have introduced notations 

G^J^\z,y) = I,iy)Ki^\z) - li-\z)K,iy), a = a,b. (22) 

In the limit 5 ^ oo the second integral on the right of (21) tends to zero (for large b/a the 
subintegrand is proportional to e^^^^/"), whereas the first one docs not depend on b. It follows 
from here that the term with the first integral in the figure braces corresponds to the Wightman 
function for the region outside a single sphere with radius a on background of the global monopole 
geometry. As a result the Wightman function in the region between two spheres is presented in 
the form 

(0|^(x)^(x')|0) = (0(«)|</,(x)¥'(x')|0(«))-^V— t^C,"/^cos0) (23) 



I 





°° zQaui (az, bz) 



dz=S^G(«) {az, Tz)G^^) [az, r' z) cosh - m?it' - t) 



where 



m \/ z'^ —W? 

QaAa^, bz) = . ) , (24) 

Ki, ' {az)li, ' [bz) - Ki, ' {bz)Ii} ' \az) 

and |0(")) is the amplitude for the vacuum state in the case of a single sphere with radius a. The 
Wightman function (0^"^ |(/7(x)(y9(a;')|0^''^) is investigated in a previous paper [12], and below we 
will mainly concentrate on the terms induced by the presence of the second sphere. Note that 
in the coincidence limit, x' = x, the second summand on the right hand side of (23) will give a 
finite result for a < r < 6, and is divergent on the boundary r = b. In Ref. [12] the Wightman 
function for a single sphere in the global monopole spacetime is presented in the form 

(0W|<^(x)^(x')|0(«)) = {Qr.MxMx')\0^) + {^{xMx'))'^\ (25) 

where jOm) is the amplitude for the vacuum state in the case of the boundary- free global 
monopole geometry. The expressions for the boundary induced part {ip{x)ip{x'))^^ for both 
regions inside and outside a single shell are given in [12]. 
Using the identity 

n^,{az,bz)G<i^Haz,rz)G^J^Hazyz)\^zl = ^^K,{rz)K,{T' z) 

Kr{az) 

-^$^Urz)lAr'z), (26) 

it can be seen that for the case of two spheres the Wightman function in the intermediate region 
can also be presented in the form 

r_/) cxD 

{OMxMx'm = (0W|^(x)<^(x')|0W)-^5;— t^Cr/^cose) (27) 



f 

Jm 







^^ z^ b^^jazM Qjh) ^^^^ ^^^Q{b) ^^^^ ^/^^ ^Qgj^ _ ^2(i' _ 

yz"^ — w? L 
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with {0^'^^\(p{x){p{x')\0^^^) being the Wightman function for the vacuum inside a single sphere 
with radius b, and 

iaz)/l!^'\bz) 



Qbi,{az, bz) 



(28) 



i?^") {az)li'^ (bz) - Rj^J'^ {bz)li''^ (az) ' 

Note that formula (27) can be also derived by the procedure described above for (23), if in 
expression (6) for the eigenfunctions we replace the function gi,^{\a, Xr) (given by (8)) by the 

function J^^ {Xr)Y^P (A6) — jj)!^ {\b)Yl^^ {^f). In the coincidence limit, the second summand on the 
right of formula (27) is finite for a < r < 6 and diverges on the boundary r = a. It follows from 
here that if we write the regularized Wightman function in the form 

{0\^{xMx')\0) = {O^MxMx')\0^) + {^{xMx'))^^^ + {ip{xMx'))'i'^ + AW{x,x'), (29) 

then in the coincidence limit the "interference" term AW{x,x') is finite for all values a < r < b. 
Using the formula for the Wightman function for a single sphere (see [12]) and equation (23), it 
can be seen that this term may be presented as 



AW{x,x') = 



where 



-a 



l-D 



TmSo 



E 



21 + n n/2 



(cos 



Jm 



zdz 



m- 



:iy("'') (r, r') cosh Wz^-m'^{t' - t) 

(30) 



r, r 



l';^iaz)K^ 



{bz) 



i^f {bz) Rlf {az) 



Glf{az,zr)Gl"'{bz,zr') 
Kl';\az)li'\bz) - K^^}{bz)iif{az) 



I,,{zr')K,,{zr) 



(31) 



Formula (29) presents the Wightman function in the region between two spheres, a < r < b. In 
the regions r < a and r > b, the Wightman functions are given by the formulae corresponding 
to a single sphere with radius a and b respectively, and are investigated in [12]. 



3 Vacuum expectation value for the field square 



The VEV of the field square is obtained computing the Wightman function in the coincidence 
limit x' — X. In the region between the spheres, from (23) and (27) we obtain two equivalent 
forms 



{0\ip'^{x)\0) = (0("V^(^)|0^°^) - 



TrnSor 



^ dz'-^k^^^G^^)\az,rz), (32) 



z" 



with a = a,b, and 



Di = {2l + D- 2) 



T{l + D-2) 



(33) 



r(D- 1)/! 

being the degeneracy of each angular mode with given I, and r(a;) is the gamma function. Using 
formula (29), the VEV for the field square can be also presented in the form 

{0\^\x)\Q) = {0^\v\x)\0^) + {v\x))^^^ + {v\x))^^^ + {v\x))^^'\ (34) 

where {<f^{x))^^ is the VEV induced by a single sphere with radius a, the "interference" term 
is given by the formula 



{^\x)) 



(ab) 



zdz 



li':\az)Ki'\bz) 



l=G 



^'''^''''■G^}{bz,rz)~ 



V^^^^Kif{az)li';\bz) - Ki';\bz)li'^\az) 

Ku,{rz) 



lS\bz) 



Kif{az) 



az, rz] 



(35) 
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and is finite for all values a < r < b. In the case of the Dirichlet scalar, by using the relation 
Iv{x)Ki,{y) > Ii,{y)Kp{x) for y > x, it can be seen that {(p^ {x))^'^^^ > 0. Note that for this 
boundary condition the both boundary induced terms {(f^{x))']^\ a = a,b, inside and outside of 
a single spherical shell are negative [12]. 

Now let us consider the limiting cases of formula (35). In the limit a ^ the subintegrand 
behaves as a^'^', and the dominant contribution to ((/^^(a;))'^"''^ comes from the summand with 
Z = 0. Using the standard formulae for the Bessel modified functions [42], noting that Ba = Ba/a, 
and assuming i^o > 0, to the leading order one has 



(ab) 



1-D 



TrSDr^-uoV^iuo) n + 2uo V2y 



where we use the notation 



x7^o(rz) 



K^\bz) 

1, 
-1, 



0, 



for 
for 



Ba = 
Ba^O 



(36) 



(37) 



For f = and a ^ 0, the dominant contribution behaves as {ip'^{x))^°'^^ ~ 1/lna. Note that 
for a minimally coupled scalar (^ = 0) one has uq = D/2 — 1. 

In the case of a massless scalar the asymptotic behavior of "interference" part (35) for large 
values b and fixed a and r can be obtained by changing the integration variable to y = zb and 
expanding the subintegrand in terms of a/6 and r/b. The dominant contribution for the sum- 
mand with a given I has an order [a/b)"^^^ (assuming that Aa / Bai^i) and the main contribution 
comes from the I = term. The leading term for the corresponding asymptotic expansion can 
be presented in the form 



2\(ab) 



2a^-^r~"- 



a \ 2^^0+1 Aa + BaUo 



TTSDaF^il^O + 1) \2bJ Aa - BaVQ 



Ag + BaVQ 
Aa - BaUo 



iT - 



f 

Jo 



dzz 



2vo 



nK^g{z) -2Sq7 zK' {z) 



nUoiz) - 2Sf^^^zi;,^{z) 



(38) 



Hence, the "interference" part of the VEV of the square of the field operator vanishes in both 
limits a — ^ and b oo. 

Now we turn to the limit a,b ^ oo for fixed b — a and a. In this limit expression (35) diverges 
and the main contribution comes from large values of I ~ o"[2(l — a/b)]^^. Introducing in (35) 
a new integration variable y = bz/vi, we can replace the Bessel modified functions by their 
uniform asymptotic expansions for large values of the order [42] . To the leading order this gives 



where we have introduced the notations 



D-2 



mi 



Cb{y)e^y/ca{y) - 1 



(39) 



mi = {b - a)^/rv? + {l/baY, Ca{y) = ^^"{^f — a = a,b, 



and 



F{y,r) 



Aa - yBa/ {b-a)' 



1 ( T — \) \ ( CL — T 

exp I 2y- ) + Ca{y) exp 2y- 1 - 2. 



Ch{y) 



b — a 



b — a 



(40) 



(41) 
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In the limit under consideration in (39) we can make the replacement 



Further, introducing instead of y a new integration variable u = ^ y'^ — t"^ — mf and converting 

to polar coordinates on the plane {u,t), the angular part of the resulting integral is easily 
evaluated. Using the standard relations for the gamma function, one receives 

(y^(x))(°'') ~ - , \n/o . , . dy 7 , ^ °\ , 7F(y,r), a, 6 ^ oo, 6 - a = const. 

^ (47r)«/2r(D/2) J^, ^ Cb{y)e''y/ca{y) - 1 ^' 

(43) 

This leading term of the corresponding asymptotic expansion docs not depend on the parameter 
a and coincides with the corresponding quantity for two parallel plates with the separation b — a, 
on background of the Minkowski spacetime. 

And finally, consider the "interference" term (y?^)^"*) in the limit of strong gravitational field, 
(J ^ 1, for fixed values a, b, r. For ^ > 0, from (10) one has S> 1, and after introducing in (35) 
a new integration variable y = z jvi , we can replace the modified Bessel function by their uniform 
asymptotic expansions for large values of the order. The integral over y can be estimated by 
using the Laplace method. The main contribution to the sum over I comes from the summand 
with I = 0, and to the leading order we receive 

^ "':°:';''r:'''':i"''^f' . v^^. <'«i. (44) 

(27ri/)^/^r"D£)VO — a ^ 

Note that, using the corresponding expressions for single sphere parts [12], we can see that in 
the limit under consideration 



and the "interference" term is suppressed compared to the single sphere contribution. For ^ = 
and a <^ 1 for the terms with / ^ one has vi ^ 1. The corresponding contribution can be 
estimated by the way similar to that in the previous case. This contribution is exponentially 
small. For the summand with Z = to the leading order over a we have ui = n/2 in (35), 
and ((^2) ~ a^"^. As we sec, the behavior of the "interference" part in the VEV of the field 
square in the limit of the strong gravitational field is essentially different for minimally and 
and non-minimally coupled scalars. This behavior is clearly seen from figure 1 where we have 
plotted the dependence of the "interference" term ((/^^(x))^"^) on the ratio r/b in the cases of 
conformally = $,d, left panel) and minimally (^ = 0, right panel) coupled massless Dirichlet 
scalars in D = 3 dimensions for a/6 = 0.5. The separate curves correspond to the values a = 1 
(a), a = 0.4 (b), a = 0.2 (c). 



4 Vacuum expectation values of the energy-momentum tensor 
and the interaction forces between the spheres 

In this section we will consider the VEVs for the energy momentum tensor operator in the region 
between two spheres on background of the global monopole spacetime. Using the standard 
classical expression, we can write the vacuum energy-momentum tensor as the coincidence limit 

{0\Tikix)\0) = lim did'k{0\^{x)ip{x')\0) + 



ik 



{0\ip\x)\0), (46) 
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0.6 0.7 0.8 0.9 




Figure 1: The "interference" part b^^^{(p'^{x))^"'''^ as a function on r/b in the cases of conformahy 
(left panel) and minimally (right panel) coupled massless D = 3 Dirichlet scalars for a/b = 0.5. 
The curves are plotted for cj = 1 (a), a = 0.4 (b), a = 0.2 (c). 



where for the point-like global monopole spacetime the nonzero components of the Ricci tensor 
are given by expressions 



9 9 ' 



(47) 



with the indices 2,3,. . . ,D corresponding to the coordinates 61,62, ... respectively. Substi- 
tuting the Wightman function (23) into (46), we obtain that the vacuum energy-momentum 
tensor has the diagonal form (as expected by the symmetry of the model) 



(OII^'^IO) = diag {e, -p, -p±, -pi_) , 



(48) 



where the vacuum energy density £ and the effective pressures in radial, p, and azimuthal, p±, 
directions are functions of the radial coordinate only. Using the Wightman function from (23) 
and the VEV for the field square from (32), the components of the vacuum energy-momentum 
tensor can be presented in the form 



q{a,b,r) =q{a,r) + qa{a,b,r), a < r < b, q = e,p,p±. 



(49) 



where q{a, r) are the corresponding functions for the vacuum outside a single sphere with radius 
a. In (49) the additional components are in the form 



qa{a,b,r) = - 



a 



l-D 



27rr"S'D 



1=0 



^ ^^z^nau^{az,bz) 



F(f rGW(az,y),G(«)(az,y)j _ ,(50) 



y=zr 



P, P±, 



where for arbitrary functions f{y) and g{y) the functions F^,'^^ [/(y))5(y)] are defined by the 
relations 



<^ [fiy),9iy)] 
f!^^ [f{y),g{y)] 

FtHfiy),9iy)] 



(1 - 40 



2y 



fg 



/V + ^(/.)'-|^i + 

(4e-l)/V-|(/5)' + 



iyf+in/2 



4^-1 + 



1-4^ 
fg, | = 4(n + l)^-n 



(n -I- l)y^ 



fg, 



(51) 
(52) 

(53) 
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and the prime denotes the differentiation with respect to y. Quantities (50) are finite for a < 
r < b and diverge at the surface r = b. Similar to the Wightman function, the components of 
the vacuum energy-momentum tensor for a single sphere case can be presented in the form 

Q{a,r) = qm{r) + qbia,r), (54) 

where qmir) are the corresponding quantities for the boundary-free monopole geometry and the 

expressions for the sphere induced parts 55(0, r) are given in [12]. 

On the basis of formula (27), the vacuum energy-momentum tensor components may be 
written in another equivalent form: 

q{a,b,r) =q{b,r) + qb{a,b,r), a < r < b, q = e,p,p±, (55) 

with q{b, r) being the corresponding components for the vacuum inside a single sphere with 
radius b. Here the additional components are given by the formula 

1^(-M = -£^12^1 dz '' 7/""'7V (f) [Gg)(6.,,),Gg)(6.,y)] ,(56) 
2Trr^bD ^ Jm V z^ - ' L ' J y=zr 

q = £,P,P±- (57) 

This expressions are finite for all a < r < b and diverge at the inner sphere surface r = a. 

It follows from the above that if we present the vacuum energy-momentum tensor components 
in the form 

q{a,b,r) = q^ir) + qbia,r) + qi,{b,r) + Aq{a,b,r), a < r < b, (58) 
then the quantities 

Aq{a, b, r) = qa{a, b, r) - q{b, r) = qb{a, b, r) - q{a, r) (59) 

are finite for all a < r < 6. Near the surface r = a it is suitable to use the first equality in (59), 
as for r — > a both summands are finite. For the same reason the second equality is suitable for 
calculations near the outer surface r = b. Using formula (30) for the corresponding part of the 
Wightman function, it can be seen that the following formula takes place for the "interference" 
parts 

A,(a6r) = -^f D, T ^ 1^ (a.)^g^ (6.) 

^'^^'^'^'^^ 2.rnSnjr:^ Vz^^ l(<:\az)Ki!;\bz) - lS\bz)K^\az) 



Fif)[/.(y),Gg)(6^,y)] F!^\K,{y),G'^\az,y)] 



I lS\bz) Kif{az) 

It can be checked that this quantities satisfy the covariant continuity equation 

dAp 



(60) 

y=zr 



dr 



+ (£>-l)(Ap- Ap^) = 0. (61) 



Note that the ambiguities of the renormalization procedure for the VEV of the energy-momentum 
tensor in the form of an arbitrary mass scale (see [18]) are contained in the boundary-free parts 
qra{f) of the corresponding components. The boundary induced parts %(«, r), a = a, 6 and 
Ag(a, b, r) are unambiguously defined for a < r < b. In particular, for the massless conformally 
coupled scalar they contain no conformal anomalies and are traceless. 
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Now wc turn to the investigation of limiting cases of formula (60). In the limit a ^ 0, the 
subintegrand behaves as o^*^' and the dominant contribution comes from the 1 = term: 



Aaia b r) « - ^^""W^fc)^" n + 2r,.uo f 



27rSDr"63i.or2(i.o) n + 2vQ Jmb V^^^^^^ I^^z) 



y=zr/b 



assuming that i^q > 0. For i/q = and a — > the "interference" parts (60) behave as 1/lna. 
Consider the limit 6 — *• oo for fixed a and r in the case of a massless scalar field. By changing 
the integration variable to y = bz and using the formula for the modified Bessel functions in the 
case of small values of the argument, we see that to the leading order the subintegrand in (60) 
behaves as (a/b)'^'^''^^ (assuming that f; > 0). The leading contribution is due to the I = term 
and one has: 



Ag a, 6, r ^ ^-^^ — — - {2vq + n)fll> - + 2f^J 

^ nK„Jz) -2SnA zK'(z) 



with notations 

/i? = Kl-40, /i? = -|, ft^ = "-^^' (64) 
/2?=0, /2? = -n-l, ft^ = 4n{n + l)Ua'. (65) 

As we see in both limits a ^ and 6 — > oo the "interference" parts for the energy-momentum 
tensor components vanish. 

In the limit a,b ^ oo for fixed values b — a and a, by the calculations similar to those for the 
"interference" part of the field square, one receives 

Aq{a, b, r) --J^-l^-^I^ dy 7 !"°r/' . AF^(t/, r), (66) 
' ' ^ (47r)^/2r(Z^/2 + 1) ^Cb{y)e-'y/ca{y)-l ^' ^ ^ 

where 

AF^{y,r) = -AF^^(y,r) = 1 - ^^^^ ~ [F(y, r) + 2], AF^' = 1, (67) 

Z[y — niQ) 

with the function F{y,r) defined by (41). These expressions are exactly the same as the cor- 
responding expressions for the geometry of two parallel plates on the Minkowski background 
investigated in [32] (note that in [32], the notations Ba/Aa = Pi and Bi,/Ai, = —(32 are used) 
for a massless scalar and in Ref. [43] for the massive case. 

And finally, let us consider the limit of the strong gravitational field, corresponding to o" <C 1. 
In this limit for ^ > one has 3> 1. After the change of variable to y = bz/ui, we replace 
in (60) the Bessel modified functions by their uniform asymptotic expansions for large values of 
the order. The leading contribution comes from the I = summand and we have the following 
limit for the "interference" parts of the vacuum energy- momentum components: 
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where v is defined in (44). In this hmit the "interference" parts are exponentially suppressed 
with respect to single spheres contributions. 

Now we turn to the interaction forces between the spheres. The vacuum force acting per 
unit surface of the sphere r = a, a = a, 6, is determined by the ^-component of the vacuum 
energy-momentum tensor at this point. By virtue of relations (49) and (55), the corresponding 
effective pressures can be presented as a sum of the pressure for a single sphere with r = a 
when the second sphere is absent and the pressure induced by the presence of the second sphere, 
Pa(o, 6, r = a). The first term is divergent due to the well-known surface divergences and needs 
additional regularization. The second term is finite and can be termed as an interaction force 
between the spheres. This additional radial vacuum pressure on the sphere with r = a, a = a, 6 
due to the existence of the second sphere can be found from (50) and (56), respectively. Using 
the relations 

G'^\rz,rz) = -Br, rz—G^::\rz,y) \y=rz= Ar, r = a,b, (69) 

oy 

they can be presented in the form 



Pa{a, b,r = a) 



1=0 



zdz 



^lauiiaz,bz) 



(70) 



[Al - 4(D - l)^AaBa - [zW + uf - nV4] 5^ } , a = a, b. 



Unlike the self-action forces, these quantities are finite for a < b and need no further regular- 
ization. For the Dirichlet scalar one has Q,ai^j{0'Z,bz) > and, hence, Paia,b,r = a) < 0. This 
means that in this case the vacuum interaction forces between the spheres are attractive. Using 
the Wronskian relation for the Bessel modified functions, it can be seen that 



[Aa - Ba{z^a^ + z/2)] ^^^{az, bz) 



^ 1 
-riaa— In 

oa 



I^^\bz) Ki^^'iaz) 



(71) 



where na = 1, = — 1. This allows us to write the expressions (70) for the interaction forces 
per unit surface in another equivalent form: 



Pa{a,b,r = a) = 



1-D 



— In 

da 



1=0 
ia) 



zdz 



n.'iaz) Kh^'ibz 



lS\bz) Kl';'{az 



A-a 



Bliz^a^ + 



a 



a,b, 



(72) 



where ^ is defined in (52). For Dirichlet and Neumann scalars the second term in the square 
brackets is zero. 

Let us consider the interaction forces between the spheres in limiting cases. For small values 
of the radius of the inner sphere, a 0, the leading contribution comes from the I = summand. 
The corresponding terms are given by formulae 



Paia,b,r = a) 
Pb{a,b,r = b) 



n + 2riaVo 



22^o-ir2(jyo)7r5i5 n + 2z/o 



dz 



z 



21/0+1 



2'^^ouor'^{uo)TTSDb^ n + 2vQ 



f 



dz- 



z 



21/0+1 



{bz) 



(73) 



[aI - A{D - l)iAbB, - [zH^ + (l/a^ - l)n(n + 1)^] B^] , (74) 



13 



where rja is defined in (37). 

For large values of the radius of the outer sphere, 6 — > oo, in the case of a massless scalar to 
the leading order we find 



Pa{a,b,r = a) 



2a^-^{Aa-BaUo) 



"2 ' a \ 21^0 



/ 

^0 



(75) 



Pbia,b,r = b) 



dz 



2vo 



a'-^{a/2bf'° Aa + BgUQ 

' TTI^or^ {uo)h''+^SD Aa - B^Uo Jo II (^) ' 



(76) 



In (76) we have assumed that Af, ^ 0. In the case = (Neumann boundary condition on the 
outer sphere) the subintegrand is equal to —{z'^ + — /^)z^'^^''[nIyQ{z)/2 — zll^{z)]~^. 

Now let us consider the limit cr ^ 1. For ^ > one has f ; 3> 1 for all I, and to the leading 
order we have 



Pa{a,b,r = a) « - 



l-D~3/2 



,-2Pln(6/o) 



2Tra^SDVb'^ - a? 



a 



a, b. 



(77) 



In the case of a minimally coupled scalar, ^ = 0, the main contribution is due to the summand 
/ = with ui = n/2 and is of an order a^~^. In this case the contributions coming from the 

I > summands are exponentially suppressed. 

As in the cases of the VEVs for the field square and energy-momentum tensor, the dependence 
of the interaction forces on the parameter a is essentially different for minimally and non- 
minimally coupled scalars. This feature is illustrated in figures 2 and 3 where we have plotted 
the dependence of the interaction forces pa{a,b,r = a) (left panels) and pi){a,b,r = b) (right 
panels) between the spheres on a/b for conformally and minimally coupled D = 3 Dirichlet 
scalars on background of the global monopole background with a = 1 (a), a = 0.4 (b), a = 0.2 
(c). 





0.1 0.2 0.3 0.4 0.5 



0.1 0.2 0.3 0.4 0.5 



Figure 2: The interaction forces between the spheres, b pa{a,b,r = a) (left panel) and 
b^~^Pb{a,b,r = b), (right panel) as functions on a/b in the case of a conformally coupled 
massless D = 3 Dirichlet scalar on background of the global monopole geometry with cr = 1 (a), 
a = 0.4 (b), a = 0.2 (c). 
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Figure 3: The interaction forces between the spheres, b pa{a,b,r = a) (left panel) and 
b^~^Pb{a, b,r = b), (right panel) as functions on a/b in the case of a minimally coupled massless 
D = 3 Dirichlet scalar on background of the global monopole geometry with a = 1 (a), cr = 0.4 
(b), a = 0.2 (c). 

5 Conclusion 

In this paper, we present the quantum vacuum effects produced by two concentric spherical 
shells in the D + 1-dimensional point-like global monopole spacetime defined by the line element 
(1). The case of a massive scalar field with general curvature coupling parameter and satisfying 
the Robin boundary conditions on the spheres is considered. To derive formulae for the vacuum 
expectation values of the square of the field operator and the energy-momentum tensor, we 
first construct the positive frequency Wightmann function. This function is also important in 
considerations of the response of a particle detector at a given state of motion through the 
vacuum under consideration [10]. The application of the generalized Abel-Plana formula to the 
mode sum over zeros of the combinations of the cylinder functions allows us to extract the part 
due to a single sphere on background of the global monopole geometry, formula (29). In this 
formula the "interference" part AW{x, x'), given by formula (30), is finite in the coincidence 
limit for all values a < r < b. The VEV of the square of the scalar field, (0|(y9^(x)|0), is given by 
the evaluation of the Wightman function in the coincidence limit. The expectation values for the 
energy-momentum tensor are obtained by applying on the corresponding Wightman function a 
certain second-order differential operator and taking the coincidence limit. In both cases the 
expectation values can be presented as a sum of boundary-free global monopole, single sphere 
induced and "interference" terms. The VEVs related to the spacetime of a point-like global 
monopole without boundaries are considered in [17, 18, 19, 20], and the effects produced by a 
single sphere are investigated in a previous paper [12]. Note that for the points not lying on 
the spheres, the boundary induced terms are unambiguously defined and the ambiguities in the 
renormalization procedure in the form of an arbitrary mass scale are contained in the boundary- 
free parts only. In particular, for the massless conformally coupled scalar the boundary induced 
vacuum energy-momentum tensor contains no conformal anomalies and is traceless. In this 
paper we concentrate on the "interference" parts in the local Casimir densities and interaction 
forces between the spheres. The application of the generalized Abel-Plana formula allows us 
to derive closed expressions for these quantities, given by formula (35) in the case of the field 
square, and by formula (60) for the energy-momentum tensor. As bulk divergences are contained 
in the part corresponding to the global monopole geometry without boundaries, and the surface 
divergences are contained in the single sphere parts, the "interference" parts are finite for all 
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values a < r < b. In particular, the integrals in the corresponding formulae are exponentially 
convergent and they are useful for numerical evaluations. We have considered various limiting 
cases of the formulae for the "interference" parts. In the limits a — > or 6 — > cxd for a fixed 
value of r, these parts vanish as a'^'^° and b~^^^°'^^^ respectively, where uq is defined by relation 
(10) with / = 0. In the limit a,b ^ oo for a fixed value b — a, the leading terms in the boundary 
produced parts do not depend on the parameter a associated with the solid angle deficit, and 
are exactly the same as the corresponding quantities for the geometry of two parallel Robin 
plates on the Minkowski background. We have also investigated the limit of strong gravitational 
fields corresponding to small values of the parameter o", a <C 1. In this limit the behaviour of 
the Casimir densities is drastically different for minimally = 0) and non-niinimally ^ 0) 
coupled scalars. In the minimal coupling case the leading terms of the corresponding asymptotic 
expansions for both field square and the energy- momentum tensor VEVs behave as a^~^ . For a 
non-mininially coupled scalar, the "interference" parts behave as {cp'^)^"'^^ ~ cj^/^~'^ exp(— 7/0") 
and Ap Ap±_ Ae/a ~ o""-^"^/^ exp(— 7/cj), with 7 = n{n + 1)^ ln(6/a). Note that in 
this case the "interference" parts are exponentially suppressed with respect to single sphere 
contributions. The vacuum forces acting on spheres contain two terms. The first ones are the 
forces acting on a single sphere then the second boundary is absent. Due to the well-known 
surface divergences in the VEVs of the energy-momentum tensor these forces are infinite and 
need an additional regularization. The another terms in the vacuum forces are finite and are 
induced by the presence of the second boundary. They correspond to the interaction forces per 
unit surface between the spheres and are determined by formula (70). For the Dirichlet scalar 
these forces are always attractive. In the limit of the strong gravitational field, a <^ 1, for 
the minimally coupled scalar field the interactions forces behave as cr^~^, whereas for a non- 
minimally coupled scalar they are exponentially small, Pa{a,b,r = a) ~ a'~^~^^'^ exp{—j/a). 
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